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We study the swelling of a flexible linear chain composed of active particles by analytical theory 
and computer simulation. Three different situations are considered: a free chain, a chain confined 
to an external harmonic trap, and a chain dragged at one end. First we consider an ideal chain with 
harmonic springs and no excluded volume between the monomers. The Rouse model of polymers 
is generalized to the case of self-propelled monomers and solved analytically. The swelling, as 
characterized by the spatial extension of the chain, scales with the monomer number defining a 
Flory exponent v which is = 1/2,0,1 in the three different situations. As a result, we find that 
activity does not change the Flory exponent but affects the prefactor of the scaling law. This 
can be quantitatively understood by mapping the system onto an equilibrium chain with a higher 
effective temperature such that the chain swells under an increase of the self-propulsion strength. 

We then use computer simulations to study the effect of self-avoidance on active polymer swelling. 

In the three different situations, the Flory exponent is now ly = 3/4,1/4,1 and again unchanged 
under self-propulsion. However, the chain extension behaves non-monotonic in the self-propulsion 
strength. 

PACS numbers: 61.25.he,82.70.Dd,61.30.Pq,87.15.A- 


I. INTRODUCTION 

In recent years, much interdisciplinary research in soft 
matter physics, fluid mechanics and biology has been de¬ 
voted to understand the motion of microswimmers in a 
low-Reynolds-number fluid M- Microswimmers can ei¬ 
ther be found as microbes such as bacteria Q, viruses 
or algae, [H or are realized artificially as self-propelled 
(“active”) colloidal particles jol-IT^. While the former 
are typically self-propelled by ch ang ing their shape, e.g. 
using beating flexible flagella [l^ , the latter are form- 
stable Janus particles exposed to a self-generated chem¬ 
ical or thermal gradient which brings the particle into 
motion The combination of self-propulsion and 

rotational Brownian motion ultimately leads to diffusion 
however with a much higher diffusion coefficient as com¬ 
pared to unpropelled (“passive”) particles 

Along the route of recent research, the shape of arti¬ 
ficial colloidal swimmers has been made more com plex 
by considering Janus-spheres [l2, l2ll-[^ . rods |2^ 

and particles of arbitrary shape |33l-[^ . At the same 
time, the flexibility of the flagella was incorporated in 
models for microbe motion j39l-l^. 

Here, we consider a flexible object that is composed of 
many active constituents, namely a linear chain of self- 
propelled particles. The motivation to do so is threefold: 
first, in a general sense, it is necessary to study how the 
collective behavior of microswimmers depends on their 
mutual coupling. Typically a pairwise interaction poten¬ 
tial is assumed between all microswimmers but what is 
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barely understood is how this behavior is affected by a 
specific strong coupling topology defining a connectiv¬ 
ity, e.g. along a linear chain. Secondly, chains of ac¬ 
tive monomers are at the interface between the physics 
of microswimmers and polymer science such that they 
establish a natural link between these two scientific dis¬ 
ciplines. It would indeed be challenging to generalize 
the broad concepts of polymer scaling theory 
towards nonequilibrium |46l-[4^ and the physics of self- 
propelled particles. Lastly, polymers of active particles 
can be realized by chaining artificial colloid^ using e.g. 
the lock-and-key technique or DNA [5l|, towards 
a linear chain of colloids. Recently, these linear col¬ 
loidal polymers have been prepared and called “polloidal 
chains” [52| or “colloidal caterpillars” [53| . When these 
colloidal beads are replaced by active colloids, the situa¬ 
tion of a flexible chain considered here is in principle ex¬ 
perimentally realizable such that theoretical predictions 
can be verified on the monomer-resolved level. Another 
realization of an active chain is a shaken granulate chain 
which has already been realized EH • Here the mil- 
limetric beads can be studied in real-space. Though the 
details of our modeling apply to microswimmers in a sol¬ 
vent, we expect qualitative similarities between an active 
polymer solution and an active granular chain. 

In previous work, appropriate models for an active 
semiflexible chain have been studied. Here a finite persis¬ 
tence length along the chain is assumed but the chain is 
typically short in the sense that it does not reach its limit 
of coiling. Some of the previous studies [56|, EH consider 
only a single active bead along the polymer, others focus 
on the dynamics of an active semi-flexible chain using 
either simulation or field theory ElHsH- Lastly, very re¬ 
cently, a one-dimensional chain of active beads (i.e. the 





case of perfect persistence) has been studied in a ratchet 
potential [6l|. 

In this paper, we focus on the swelling behavior of a 
flexible chain composed of active beads using both ana¬ 
lytical theory and computer simulation. In doing so we 
discriminate between an ideal chain (i.e. a chain without 
self-avoidance) and a self-avoiding chain. For the former 
case, we generalize the traditional Rouse model 13163 
of polymer dynamics towards the situation of active 
monomers, for the latter we use computer simulations. In 
the spirit of the simplicity of the Rouse model which as¬ 
sumes a linear bead-spring chain with harmonic coupling 
between nearest neighbors and neglected hydrodynamic 
interactions we introduce an activity for the monomers 
but neglect any correlations in the activity of neighbor¬ 
ing beads. This generalized Rouse model with no explicit 
aligning interactions for active polymers is then solved 
analytically in two dimensions. As a basic result, the 
Flory random-walk exponent z/ = 1/2 which measures 
how the typical extension of the chain scales with the 
monomer number is not affected by activity but the pref¬ 
actor is corresponding to a higher effective temperature 
when mapped onto a corresponding passive chain. This 
implies that the extension of an active chain increases 
with the strength of the self-propulsion (or Peclet num¬ 
ber). 

Self-avoidance in two dimensions yields a Flory 
swelling exponent v = 3/4 for passive monomers. Our 
computer simulations show that for long chains this ex¬ 
ponent is not affected by self-propulsion of the monomers. 
This is in line with experiments using a shaken granular 
chain Ell where the Flory exponent was found to be 
unaffected by the activity. However, as a function of the 
self-propulsion strength, there is a non-monotonicity in 
the chain extension which is absent for an ideal chain. 
Interestingly, this effect was found recently in the reverse 
set-up of a passive chain in an active bath (bH, [b^ which 
can be realized for granulates [bl|, and which shows the 
same Flory scaling for large monomer number. 

We also solve the Rouse model of an active ideal chain 
either confined to an external harmonic trap or dragged 
at one end. In these situations the Flory exponent is 
z/ = 0 and z^ = 1, respectively, and also not changed 
by activity. But again the prefactor is affected by self¬ 
propulsion. Self-avoidance leads to the Flory exponents 
z/ = 1/4 and z/ = 1 in these two situations. Our simula¬ 
tions show that these exponents are unaffected by inter¬ 
nal activity of the chain. 

The paper is organized as follows: we describe our 
model and the different situations in Sec.HIland the simu¬ 
lation in Sec. imi Results are discussed in Sec. lEl and we 
conclude in Sec. m where we comment on how an active 
polymer can be realized experimentally. 



FIG. 1. Sketch of an active polymer consisting of N monomers 
which are self-propelled with an effective driving force F along 
their orientation Un for the three studied situations: (a) free 
chain, (b) chain in harmonic confinement I/ext and (c) chain 
with a dragging force Fd applied to the first monomer. 


II. ROUSE MODEL FOR AN IDEAL CHAIN OF 
ACTIVE PARTICLES 

In our Rouse-like model we describe the polymer as 
a linear chain of N beads. Neighbors are connected by 
harmonic springs [bb| with a spring constant k. We fo¬ 
cus on the two dimensional case. At a given time t, the 
monomers are at positions rn{t) = [xn{t),yn{t)] where 
n = 1,...,V, and possess orientations described by the 
unit vectors Un = (cos^n, sin^^). The self-propulsion 
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velocity of each monomer is introduced via an effec¬ 
tive driving force acting along this orien¬ 

tation, see Fig. [IJa). A single monomer experiencing 
overdamped dynamics with a translational friction coeffi¬ 
cient 7 would then move with the constant self-propulsion 
speed vq = F/^ We neglect any hydrodynamic in¬ 
teractions as in the traditional Rouse model [62|. The 
overdamped equations of motion for the positions rn{t) 
of the n-th bead of the chain can then be written as 


-‘S'« 

whereby d‘^rnldn^ = (r^+i + Vn-i — 2 rn) corresponds 
to the spatial differential quotient of neighboring beads, 
using the constraints ro = 0 and rw+i = 0. Assuming 
the Stokes-Einstein relation applies not too far from equi¬ 
librium, the translational friction coefficient 7 is close to 
k^T!D where D is the short-time diffusion coefficient for 
a single bead and /cbT the thermal energy. The random 
forces ^n{t) are Gaussian distributed with zero mean and 
variance (|(t) (g) = 26{t — t')l (/cbT)^ /D. The ori¬ 

entation of the particle is described by the rotational 
Langevin equation 


dUn _ . . . 

7r ~ Cn ^ • 


( 2 ) 


Here ^ is a Gaussian-distributed torque with zero mean 
and variance {C{t) G — t')l (/cbT)^ /Dr and 

7 ^ is the rotational friction coefficient, which is close to 
k^T/Dr. The rotational diffusion coefficient Dr is deter¬ 
mined from the relation D/Dr = 4:b‘^/3. 

In our model, b = ^Jk^T/k characterizes a typical 
bead extension or bead distance generated by the finite 
temperature, the so-called Kuhn length. In terms of this 
length, a typical elastic energy contained in the harmonic 
spring is given by kb^ such that we define the ratio be¬ 
tween thermal and elastic energy as the parameter 


3k^T ' 


(3) 


The strength of self-propulsion, on the other hand, is 
characterized by the Peclet number 



(4) 


We do not consider explicit aligning interactions between 
neighboring beads. 

In this paper we focus subsequently on three different 
situations, namely a free chain, a chain confined to an 
external harmonic trap, and a chain dragged at one end. 
These situations are sketched and summarized in Fig. [TJ 
A symmetric external harmonic trap potential is shown 
in Fig. [IJb). It is given by 


^ext ('^) — 2 ’ 


( 5 ) 


and can be incorporated by adding the corresponding 
external force —KTn to the right-hand-side of Eq. O- 
This introduces a further dimensionless parameter into 
the model, namely the ratio of the external and internal 
spring constants njk. 

Lastly, a constant drag force F^ — —F^ex is applied 
to the first monomer of the chain, see Fig. me). The 
dragging force results in an anisotropic mean shape of 
the chain. The relative strength of the dragging can be 
described in terms of elasticity determined by the dimen¬ 
sionless parameter Fd/kb. 


III. COMPUTER SIMULATIONS FOR A 
SELF-AVOIDING CHAIN 


In the Brownian dynamics simulations, we model the 
active polymer as a sequence of N coarse-grained spring 
beads in analogy to previous works considering short flex¬ 
ible rods [gsI, l69|. For simplicity, interactions between 
the active monomers are modeled by a smooth repulsive 
WCA (Weeks-Chandler-Andersen) potential 


UwcA{r) = 




r < 

r > 2 ^/^. 


(6) 


Here a denotes the diameter of a single bead and e = /cbT 
is the interaction strength. These quantities represent 
the length and energy units, while times are conveniently 
measured in units of the Brownian time r = a^/D. The 
self-avoidance of the chain is incorporated via the finite 
length a. 

Springs between neighboring beads are introduced via 
a so-called FENE (finitely extensible nonlinear elastic) 
potential 0 


^FENE(Gi) = 



(7) 


with neighboring beads i^j and their distance Vij = |r^ — 
Tj\. The spring constant is fixed to K = and the 

maximum allowed bond-length to Rq = 1.5cr. 

In our Brownian dynamics simulations, the same over¬ 
damped equations of motion [see Eqs. o and @] were 
integrated with a small finite time step At, but now with 
self-avoidance and FENE-chains. The external harmonic 
potential and the dragged monomer were included us¬ 
ing the same expressions as introduced in the previous 
chapter. 

Statistics are gathered for up to 50 independent initial 
configurations along times of t = lO^r after an equilibra¬ 
tion period. The time step used in our Brownian dynam¬ 
ics simulations is At = 10 “^r. All our simulations are in 
two spatial dimensions. 
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IV. RESULTS 
A. Free ideal chain 

First we focus on a free ideal chain which is long {N ^ 
1). In our analytical calculations, we consider the linear 
transformation 


/ N 

djTi (f)pjiVji{t^ ( 8 ) 

and choose the coefficients (j)pn such that the equation of 
motion for Xp has the form 

lp^ = -KXp + Fp{t). (9) 

This leads to the relations 

-Fp = — / drKppnFn , (12) 

7 Ji 

using Fn(t) = {Fiin + Cn(^)) • Ca and jp = 2N^, with 
p = 1,2,... and a = x,y. These equations can now 
be applied to calculate the mean square of the gyration 
radius as 

dn{rn-rcm)a'^ 

oo 

= 2Y,{Xp{t)Xp{t)). (13) 

p=l 

Here the center of mass is given by Vcm = ^ Yln=i 
and (• • •) denotes a noise and time average. The end-to- 
end distance is given by 




0 5 10 15 20 25 30 

(b) Pe 


FIG. 2. (a) Reduced end-to-end distance Re (solid lines) and 
radius of gyration Rg (dashed lines) for a free active polymer 
chain at Pe = 10. (b) Relative change in induced by 

activity Pe for fixed number of active beads N. 


Interestingly, the Flory exponent describing the basic 
scaling of spatial extension of the active polymer, see 
Eqs. and (pT|) . is the same as for a passive chain, i.e. 
we obtain for large N both 

{R%) ~ N'^'^ (18) 


REocit) = rNa(t) - riait) = “4 ^ 4 ) 

p: odd int 

for each component, where the sum is over all odd values 
of p. 

For a long chain (V 1) of self-propelled particles we 
therefore find 


{Xp{t}xs}) 

{Rl)/h^ 


Po 


7p7 


2 kp “h D-p'^p kp 




2 N 
9 T 
AN 

sT 


2 o 

l + -Pe^ 

2 o 

H--Pe^ 




(16) 


(17) 


and 


{R%) ~ N'^'' (19) 

with z/ = 1/2. This implies that for long chains with¬ 
out aligning interactions the effect of self-propulsion is 
not so dramatic that the scaling is changed. We expect 
this to hold true even for short-ranged aligning interac¬ 
tions. Long-ranged aligning interactions, however, could 
have an influence on the Flory scaling exponent, similar 
to the scaling of polyelectrolyte chains with long-ranged 
Coulomb interactions 0 . 

More precisely, the expressions m and dm) are iden¬ 
tical to those for a passive polymer at a higher effective 
temperature 


Teff/T = 1 + |pe2 . 


( 20 ) 
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The same mapping onto an effective temperature has 
been made for a single self-propelled monomer in a grav¬ 
itational field 0 and has been tested for a single self- 
propelled bead along a chain [56|, • 

Finally, the long time diffusion coefficient for the chain 
is given by 


An analytical calculation for the ideal chain is still pos¬ 
sible. However, the external potential leads to a shift in 
the values of kp^ which are now given by 



Dl 


D 

N 


l + fpe 



( 21 ) 


Accordingly, the radius of gyration and the end-to-end 
distance are now given by 


which is identical to the case of a single active colloid 
but reduced by a factor 1 /A^. 

We now turn to effects of self-avoidance. In two spa¬ 
tial dimensions self-avoidance changes the Flory expo¬ 
nent from 1/2 to 3/4. Let us briefly recapitulate Flory’s 
argument: The free energy of a passive, self-avoiding 
chain in two dimensions is composed of two parts 

E ^ F’ex + Eq\ . (22) 

The first one, F’ex arises due to excluded 

volume effects and leads to chain swelling. The second 
term Eq\ ^ kR^ /N incorporates the elastic properties of 
the chain and leads to chain shrinkage. In equilib¬ 
rium the total energy becomes minimal, leading to the 
well known Flory exponent v = 3/4. 

We confirm by simulations that the scaling exponent 
of the spatial extension, measured by {R%) or {R%), is 
indeed 3/4 unaffected by the activity, see Fig. [2fa) |^ . 
This is in agreement with the experimental findings of an 
active granular chain [HJ, where for chain lengths N < 
128 the Flory exponent of a self-avoiding chain = 3/4 
has been found. 

A striking difference generated by self-avoidance, how¬ 
ever, is the Peclet number dependence of the chain ex¬ 
tension. As documented in Fig. Hb), we find a non¬ 
monotonicity of the chain extension for increasing self¬ 
propulsion. The chain first shrinks and then swells. 
The initial chain shrinkage can be attributed to self¬ 
propulsion of a particle encaged by its hard neighbors. 
Activity will on average drive the particle away from its 
neighbors which will contribute to a zig-zag configuration 
of the chain implying a shrinkage. Increasing the activity 
further will finally swell the chain since the larger driv¬ 
ing force will extend the chain. The minimum in the 
chain size therefore occurs when the self-propulsion force 
E is getting comparable to the elastic force acting on a 
particle. Interestingly, very recently, this effect has ben 
found in the reverse set-up of a passive chain in an active 
bath [gsl, Q showing that an active bath particle can be 
viewed as forming a joint unit with the polymer chain 
once it is colliding with the chain. 


B. Harmonically confined chain 

Now, we include a confining harmonic potential, see 
Eq. ([5]). We vary the prefactor k. with respect to the 
spring constants keeping the activity of the monomers 
fixed. 




coth N 


1 

TV 


Xn 


Xu E k 

2 k 
3A K, 


Xhv E k 

Xk 


coth N 


Xk E k 


Xk 




coth N 


1 

N 


(24) 


and 


{Rl)E = ^Pe^ 


Xk IN XK + k 







( 25 ) 


respectively. In the limit of large N we find the Flory ex¬ 
ponent z/ = 0 which implies that the chain is completely 
localized. Again, the activity of the monomers only af¬ 
fects the prefactor but not the scaling behavior of the 
whole chain. As a function of monomer number AT, a 
plateau in the spatial extension can be observed above a 
threshold chain length Nf which scales with the strength 
of the applied external potential as 

"• = + s’’”’) ' ■ <“> 

Consequently, for large self-propulsion strengths, Pe 1 , 
the threshold chain length decays as Nt ^ Pe~^. Below 
this monomer number the well known Flory exponent 
u = 1/2 is found, see Fig . [3l^a), reproducing the trends 
for a passive polymer |75| . 

For a self-avoiding passive polymer chain, Flory’s ar¬ 
gument, see Eq. (|2^ . can be extended by including an 
external harmonic potential to the total free energy as 

E ^ Eex + Eel + Eext 5 (27) 

where the scaling of Eext can be obtained by integrating 
a constant sharp-kink monomer density profile inside a 
circle of radius R as Eext ^ kNR^. The latter two con¬ 
tributions will enforce a compactification of the polymer 
chain. Eor long chains (AT ^ 1 ) the external energy Eext 
dominates the elastic energy which leads to a new expo¬ 
nent u = 1/4. Below a threshold monomer number Nt 
the elastic properties of the chain dominate such that the 
Elory exponent is z/ = 3/4. Again this monomer number 
scales as Nt r\j ^Jkfn. 
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C. Dragged chain 


We will now treat the situation of a constant drag force 
Fd applied to the first monomer of the chain, see again 
Fig. IDJc). The equations of motion are the same as in 
Eqs. m and 0 but now include the term —FdCx of 
the right hand side of Eq. o for The analytical 

solution of the Rouse model can be extended to an ideal 
dragged chain. In detail, we obtain 




1 

y 2 kp T Dx'^p 

4 E/ 

H- — , 

kpkq 


kp kp J 


^pq 


p,q>l. (28) 



FIG. 3. Reduced end-to-end distance Re (solid lines) and 
radius of gyration Rg (dashed lines) for a harmonically con¬ 
fined polymer in the case of (a) no self-avoidance and (b) a 
self-avoiding chain for fixed activity Pe = 10 of the monomers 
and varied confinement strengths k. (c) Threshold number of 
monomers Nt as a function of n/k. 


This implies that 

{«l>/^.^=5yl+JPeyi(§)^^ (30) 


such that for long chains {N ^ 1) the Elory exponent is 
z/ = 1 in the dragged situation, see Eig. HJ^a). Intuitively 
this is expected as a strong drag force stretches the chain 
along its pulling direction Cx . As becomes directly visible 
from the expressions (|29|) and PH) . there is a crossover 
in the scaling of the squared chain extension from N to 
N‘^. In fact, for small drag forces and small N the Elory 
exponent is still that of an undragged chain (z/ = 1/2). 
Eor large self-propulsion (Pe ^1), the crossover between 
the two regimes in N occurs at a threshold of 




rPekb\ 

) 


( 31 ) 


Our computer simulations for a self-avoiding chain con¬ 
firm the same qualitative behavior. Now there is a 
crossover from the Elory exponent z/ = 3/4 of an ideal 
self-avoiding chain to the stretched case z/ = 1, see 
Eig. m^b), whereby the crossover number of monomers 
Nc decreases as given in Eq. m, see Eig. m^c). Again, 
the activity of the chain does not affect the scaling ex¬ 
ponent which stays to be z^ = 1 as known from a passive 

chain [zi[z3- 

Einally, in Eig. O we study the Peclet number depen¬ 
dence of the extension of a dragged chain. Indeed the 
threshold chain length for the crossover Nc scales as pre¬ 
dicted in Eq. m and shifts for higher self-propulsion to 
higher N, see Eig. [Sfc). 


In case of a self-avoiding active polymer chain, we find 
the same Elory exponent z/ = 1/4 in the limit of large 
N, see Eig. [3fb). As before, the scaling behavior is not 
drastically altered by the activity. Eor varied confinement 
strenghts we can confirm the predicted scaling for the 
threshold monomer number Nt for an ideal as well as a 
self-avoiding chain of active particles, see Eig. He). 


V. DISCUSSION AND CONCLUSIONS 

We have generalized the traditional Rouse model of 
polymer dynamics to the situation of an active polymer 
consisting of self-propelled monomers. By analytically 
solving the Rouse model and numerical simulations of a 
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(a) N 




(c) Fd/{kb), Fj{Kcj) 

FIG. 4. Reduced end-to-end distance Re (solid line) and ra¬ 
dius of gyration Rq (dashed line) for an active polymer which 
is dragged with a force Fa in the case of (a) no self-avoidance 
and (b) a self-avoiding chain for a constant activity Pe = 10. 
(c) Crossover number of monomers Nc for fixed activity as a 
function of reduced drag forces Fd/{kh) respectively Fd/{Ka). 


self-avoiding polymer we have shown that the well known 
Flory scaling exponents are still valid in the presence of 
activity. For an ideal chain of active particles, the ac¬ 
tivity only affects the prefactor and can be considered 
as an effective temperature when mapped onto a corre¬ 
sponding passive polymer. In the case of a self-avoiding 
chain, the spatial extension reveals a non-monotonicity as 
a function of activity which is absent for an ideal chain. 



(a) N 



FIG. 5. Reduced end-to-end distance Re (solid line) and 
radius of gyration Rg (dashed line) for a chain of active par¬ 
ticles which is dragged with a fixed force Fd in the case of: 
(a) Fd/(kh) = 1 with no self-avoidance and (b) Fdl^Kcr) = 1 
for a self-avoiding chain for varied activity Pe. (c) Grossover 
number of monomers Nc for a constant drag as a function of 
activity Pe. 


The cases of harmonically confined and dragged active 
polymer chains reveal a crossover in scaling for the chain 
swelling but again activity does affect the Flory exponent 
for long chains. 

Chains of active particles can be prepared in nature 
by chaining artificial colloids. This can be achieved us¬ 
ing the lock-and-key technique or dipolar or patchy 














8 


colloids th at p refer chaining, for some recent realizations, 
see Refs. [^, In principle, one can use active col¬ 
loids as entities to obtain a one-dimensional chain of ac¬ 
tive particles. This has not yet been realized but is at 
least conceivable. The main advantage of the colloidal 
realization is that the motion of the chain can directly 
be observed in real-space such that the statistical aver¬ 
age needed for the mean size is obvious. Moreover the 
activity can be controlled from outside such that a di¬ 
rect comparison with a passive thermal chain is possible, 
which enables a test of the effective temperature scaling. 
Though our model is designed for Brownian motion in a 
solvent, we expect that the key trends are the same as 
for driven granulates. Therefore, it would be interesting 
to perform more granulate experiments in confinement 
and under drag to test our further predictions. A har¬ 
monic confining potential can easily be realized by shak¬ 
ing the granular not on the plane but within a paraboloid. 
Finally, dragging the first monomer can be realized by 
charging one end-monomer using electrofriction and 
placing the whole granulate chain into a homogeneous 
electric field along the vibrating substrate. 

Future research should consider three spatial dimen¬ 
sions where the description of the orientational depen¬ 
dence of Brownian motion is more complicated [33| , but 
similar results are expected. While the scaling should still 


correspond to that of a passive chain, the prefactor will be 
affected by the activity. In particular, in three spatial di¬ 
mensions the activity is less important, which is already 
visible for the effective temperature mapping. For ac¬ 
tive spheres in three dimensions Teff/T = 1 + holds 
with a prefactor 2/9, instead of Teff/T = 1 -h |Pe^ in two 
dimensions where the prefactor 2/3 is stronger . An¬ 
other line of future research concerns a chain with both 
active and passive monomers. Here, again, we expect 
the same scaling exponent but a prefactor in the effective 
temperature mapping which is smaller than in the case 
considered here, where any monomer is active. More¬ 
over, the influence of hydrodynamic interactions needs 
to be considered in more detail following the work of 
Ref. [iH which applies to short and stiff filaments. If 
the solvent is viscoelastic, new effects come into play, see 
e.g. Refs. (^,[1^, which also need to be explored in the 
future. 
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